Unification of the standard and gradient theories of phase transition 



B. I. Lev and A. G. Zagorodny 
ov Institute for Theoretical Physics, NAS Ukraine, Metrolohichna H~b, Kyiv 03680, Ukraine 

(Dated: March 8, 2013) 

We show, that the standard model of phase transition can be unified with the gradient model 
of phase transitions using the description in terms of the gradient of order parameter. The gener- 
alization of the gradient theory of phase transitions with regard to the fourth power of the order 
parameter and its gradient is proposed. Such generalization makes it possible to described wide 
class of phase transitions within a unified approach. In particular it is consistent with the nonlinear 
models, that can be used to describe a phase transition with the formation of spatially inhomoge- 
neous distribution of the order parameter. Typical examples of such structures (with or without 
defects) are considered. We show that formation of spatially inhomogeneous distributions of the 
order parameter in the course of a phase transitions is a characteristic feature of many nonlinear 
models of phase transitions. 
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Withing the context of the general theory of phase 
transitions a system treated as a continuous medium is 
assumed to have a ground state which can always be de- 
scribed in terms of the order parameter. Such order pa- 
rameter can have various geometrical presentations, for 
example, a scalar field in the case of condensed matter 
1 1, a fundamental scalar field in the quantum field theory 
2 1 , a magnetization vector in the theory of magnetism [3| , 
a second-rank tensor in the liquid crystal theory Q,etc. 

To introduce the order parameter that determines a 
stable state of condensed matter, we have to consider pos- 
sible deformations of the distribution of the field, in par- 
ticular, the disordered configuration of the ground state. 

The phase transition associated with the system with 
broken continuous symmetry can be described in terms 
of the relevant order parameter. In particular, accord- 
ing to the Landau theory the free energy density can be 
presented in terms of the order parameter as given by 



/=(V^(r)) 2 +%(r)) 



(1) 



where (f(r) is the order parameter and W(ip(r)) describes 
the order parameter dependence of the free energy that 
is assumed to be known. 

In the well-known standard model of phase transitions 



W{ V {v))= l -a^{v) + \b V \v). 



(2) 



With the dimensionless variable tp 2 (r) = h<f ' ^ being in- 
troduced the standard dimensionless free energy reduced 
to the form given by 



/ = / 2 (V^(r)) 2 + ((^(r)) 2 -l) 5 



(3) 



with the potential being written in the standard form 
and I 2 = - being the characteristic length. Making use of 
this expression for the free energy, we can find the spatial 
distribution of the order parameter and thus, describe the 
properties of the new states which can be formed after 
the phase transition. 



It should be noted that in the case of system with the 
gradient of the order parameter, it looks reasonable to in- 
troduce into the free energy functional the term responsi- 
ble for possible interaction between the order parameter 
and its gradient. This coupling can regularize possible 
perturbations of the order parameter and thus confine 
the spatially inhomogeneous state of the system. 

On the other hand, singular perturbation models, in- 
volving higher order of the the spatial derivations have 
provided a new insight on the role of additional physical 
features of the system under consideration, on the details 
of phase transitions, in particular on the way how to de- 
scribe the contribution of the surface energy. In this case 
restricting the accuracy of the free energy functional by 
the first order derivatives of the order parameter leads 
to the solution with the homogeneous distribution of the 
order parameter only. 

The main idea of the present contribution is to gen- 
eralize the phenomenological theory of phase transitions 
by introducing the second-order spatial derivatives into 
the order parameter dependence of the free energy. Such 
generalization is reasonable in view both mathematical 
and physical arguments. From the mathematical point 
of view, our argument is that the order parameter can 
be treated as a vector quantity and thus it is possible to 
rewrite the free energy (1) in the form 



/ = Z 2 (V^(r)) 2 + ((|^(r)|) 2 -l) 5 



(4) 



where the order parameter is a vector function. Such 
problem arises if we assume that the order parameter 
tp(r can be presented is in term of a gradient of some 
other scalar function Vu(r. This assumption leads to the 
knownpresentation of the free energy in the Aviles-Giga 
form @,@,i.e. 



f = l' 



(|VVu(r)|) 2 + ((|V U (r|) 2 -l) 2 



(5) 



As is known, this presentation has various physical ap- 
plications, e.g. the description of smectic liquid crystals 
thin film blisters and convective pattern for- 

mation 9]. Physically, such model can be regarded as 
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the Landau model applied to a system with vector or- 
der parameter. Some well-known Landau theories have 
similar features. For example, the energy of a smectic- 
A liquid crystal has been described within such model 
[f| where represents the director field of the liq- 

uid crystal. The observed focal-conic defect structures 
can also be described in terms of this functional [6[ . The 
micro-magnetics give one more example of the applica- 
tion of the functional (5). In particular, it can be used 
to describe magninization constrained by \m\ = 1 within 
and m = outside the micro-magnetics [111 ], The un- 
known u(r) is purely curl-free, while \m\ just prefers to 
be divergence-free. Moreover, \m\ is restricted to unit 
vectors, while u(r) prefers to have unit magnitude. But 
the similarity should be clear, particularly for an isotropic 
ferromagnet 

Another motivation to use the functional (4) arises 
from recent phenomcnological modeling of blisters in 
compressed thin 1ms [7j. Early it has been suggested 
that the fold patterns of such blisters could be described 
by minimizing the sum of membrane and bending ener- 
gies With some simplication, this problem can be 
reduced to the free energy, given by (4). The same free 
energy can also be obtained for an equilibrium state of a 
free surface. In this case u(r) represents the height pro- 
file of the sheet (relative to a flat reference state). This 
model describes a fluctuating fluid membrane. 

The interpretation of the free energy appears also in 
the phase diffusion theory of pattern formation proposed 
by Cross and Newell |14j],[l5l- Of course, the physical 
justification of the models discussed above can be crit- 
icized [13] ■ Nevertheless, they demonstrate the need to 
introduce new representations for the free energy func- 
tional in term of the order parameter. We have referred 
informally to the existence of an asymptotic variational 
problems. Now, let us consider the possibility to general- 
ize the representation of the free energy to the case when 
it depends not only on the order parameter, and its gra- 
dient, but also on their combination. Such generalization 
looks quite reasonable if one bears in mind that, in the 
case of the functional with the first order spatial deriva- 
tive, the order parameter describes a phase transition in a 
spatially homogeneous system, while the presence of the 
second-order derivative makes it possible to describe the 
formation of spatially inhomogeneous structures. Thus, 
we can expect that the combination of the gradient terms 
with the scalar order parameter can be responsible for 
the self-consistent influence of the order parameter on 
the spatial distribution and parameters of the ordered 
structures. In other words, the coupling between the or- 
der parameter and its derivatives can influence possible 
deformations and confine possible inhomogeneous stable 
structures of the system. 

So, let us postulate, that the order-parameter depen- 



dence of the free energy is given by 

/ = aZ 2 (|V^(r)|) 2 + 6((|^(r)|) 2 - 1) 2 + 
+ml 4 (|VV<^(r)|) 2 + n{l 2 (|V<^(r)|) 2 - I) 2 + (6) 
+ap 2 {r) {lVp{r)f 

where I is the length of the order-parameter changing, 
a,b,m,n and c are parameters describing the influence of 
the gradient order parameter and coupling between the 
order parameter and the appropriate derivative of this 
parameter. Hiving introduced the operator D — ZV we 
can rewrite the free energy density as 

/ = a (|^(r)|) 2 + 6((|^(r)|) 2 -I) 2 + 
+m (\D 2 tp(r)\) 2 + n(l 2 (\D V (r)\ f - 1) 2 + (7) 
+ap 2 {v){Dp{v)f 

If all the coefficients except a and b are equal to zero, 
we come to the standard theory of phase transitions. 
If a, b and c vanish wee come to the standard gradient 
theory. In the case of large values of m ( m is larger 
than other coefficients) we obtain the eikonal equation 
V<y9(r) — I which we supplement with the boundary con- 
ditions ip(r) — at the boundary. This eikonal equation 
has no smooth solution, but it has innitely many Lips- 
chitz solutions. We can suggest that the energy can be 
concentrated at the discontinuities of <p(r). Thus, the sin- 
gular part of the order parameter can provide a selection 
mechanism for the perturbations of the eiconal equation 
solutions which minimize the free energy. 

In the general case, the minimum of the free energy 
satisfies the Euler-Lagrange equation: 

V df df 2 df 

8<p(r) dtp(r) dDip{v) dD 2 ip{v) K ' 

which is reduced in our case to: 

mD 4 tp(r) - (a - 2n + 6n(Dip(r)) 2 + cip 2 {r))D 2 ip(r) 
-ctp{r)(Dtp{r)) 2 - 26<p(r)(l - ip 2 {r)) = 

(9) 

Let us consider in the one dimensional case some proba- 
ble solutions of the Euler-Lagrange equation for various 
combinations of the coefficients introduced. 
a) Linear solutions 

i) We shall look for a solution similar to the solution 
of the equation Dip — —<p and D 2 (p = (p. Substituting 
this solution in the Euler-Lagrange equation yields 

m(p-(a-2n+ctp 2 + 6n(p 2 )tp-c(p 3 ~2bp(l-(p 2 ) = (10) 

This equation leads us to the following relations between 
the coefficients, m — a = 2(6 — n) and b — c = 3m. In the 
standard model of phase transitions a = b = I and thus 
we have m + 2n = 3 and c = I — 3n. If m = n = 1 the 
coupling constant c = — 2. Ifa = I and b = — 1 and m = 
1 and n = — 1 the coupling constant c = 2. In the case 
of the gradient presentation a = b = the exponential 
solution can be realized only for n = 1, m = — 2 and 
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c = 3. In the case of standard model with n = m = 
the exponential solution exists for 6 = 1, a = — 2 and 
c = 1. Thus, the exponential solution is realized for 
various combination of the introduced coefficients. 

ii) Now let us look for another possible linear so- 
lution in the form of a periodical function, namely 
<p(r) = <pexp(ikr). This solution satisfies the equation 
D 2 p(r) = — k 2 p(r). Substituting of this solution into 
Euler-Lagrange equations leads to the relation 

mk A p(r) + (a — 2n + cp 2 + 6nk 2 p 2 )k 2 p(r) , . 
-cp 2 <p(r) - ^bp(r)(l - p 2 ) = (Uj 

which generates two conditions i.e. ink 4 + (a — 2n)k 2 — 
2b and 6nk 2 + 2b = 0, for which the spatial periodical 
distribution of the order parameter does not depend on 
the coupling constant. The conditions thus obtained can 
be used to find the wavelength of the periodic distribution 
of the order parameter. The second condition yields k 2 = 
— On the other hand from the first condition one 
obtains k 2 = — a ~] m n ■ This means that these conditions 
are consistent with each other if mb — 3n(a + An). If 
k = 1, when the period of the structure is equal to 
the size of the system under consideration, we obtain 
m = — (a + An) or b = — 3n. In the case of the standard 
theory of phase transitions with a = b = l,we obtain 
relation for the other coefficients i.e. n = — | and m = 
■|. It is important, that in the case under consideration 
the solution for the order parameter is independent of 
the coupling constant. The next step is to consider the 
realization of possible nonlinear solutions which satisfy 
the minimum of the free energy density. 
a) Nonlinear solutions 

i) We start from the nonlinear solution of the standard 
model of phase transitions, namely, from the well-known 
the soliton solution which satisfies the equation 

D 2 ip = 2p(p 2 - 1) (12) 

The solution of this equation is ip = tanh(lr). This so- 
lution satisfies an equation with the high-order deriva- 
tives that is given by Dp = (1 — p 2 ) and next relation 
D 4 p = 8ip(l — p> 2 )(3<p 2 — 1) Having substituted these re- 
lations into Euler-Lagrange equation one finds that thus 
solution can be realized under the conditions n = 0, 
16m + 2(a — b) — c = and c = 8m. In the case of 
the standard theory of phase transition a = b all the 
other coefficients should be equal to zero. In the case of 
the generalized model we can observe the soliton solution 
for a = 1,6 = 2, c = 2 and m — j. We can also consider 
other nonlinear solutions of the standard model. 

ii) The second example of the nonlinear solution of the 
standard model concerns the one p — sech(lr) that sat- 



isfied the nonlinear equation 

D 2 tp = p(l - 2p 2 ) (13) 

Assuming also that the equations Dp 2 = p 2 (l — p 2 ) and 
D 4 p = p(l — 20p 2 + 2Ap 4 ) are satisfies one obtains an 
equation for the coefficients of the generalized model. In 
particular, the solution in the proposed form exists, if 
m = a + 26, c = —8m and a — c + 6 = 10m. This means, 
that in the case of the standard model a — b does not 
lead to the minimization of the generalized free energy 
functional. 

iii)Finally, let us consider whether the solution p = 
lnch(lr) can be realized. This solution satisfies the equa- 
tions D 2 p = (1 - (Dp 2 ), Dp — tanh(lr) and D A p — 
2(1 - (Dp) 2 )(3(Dp) 2 - 1). It follows from the Euler- 
Lagrange equation that this solution minimizes the free 
energy functional under the conditions a = b = c = 0, 
i.e., we obtain the well- known form 

/ = (D 2 p) 2 + (1 - (Dp) 2 ) (14) 

The absolute minimum of the free energy is achieved for 
p = 1 and Dp = 1, that cannot be satisfied. This means, 
that the the eikonal solution satisfies the condition of the 
free energy divergence. Such conditions can be realized 
only for spatially inhomogeneous order parameters. 

Thus, we have proposed a generalization of the stan- 
dard and gradient theories of the phase transitions by 
introducing the coupling between the order parameter 
ant its gradient. Such generalization can be employed 
to describe the phase transitions from spatially homo- 
geneous to inhomogeneous states. It is shown, that the 
solution of the standard and gradient models of phase 
transitions can be inconsistent with the Euler-Lagrange 
equation generated by the generalized functional of the 
free energy. In the general case the requirement of con- 
sistency of the known solutions with the generalized de- 
scription can be archived by the appropriate choice of the 
coupling between the order parameter and its gradient. 

Specific examples considered in the present contribu- 
tion shows, that the structure formation observed exper- 
imentally can be described by various phenomenological 
free energy functional which correspond to various sets 
of coefficients, i.e. there is no unique functional repre- 
sentation of the free energy related to the chosen spa- 
tially inhomogeneous configuration of the order param- 
eter. This uncertainty is generated by the phenomeno- 
logical description and could be eliminated in the micro- 
scopic calculations. 
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